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INTRODUCTION 


The purpose of this stuay was to select fnonthly weather vanaDies that 
cou1o De usea to estimate yields of tne Aryentina soybean crop. The soybean- 
growing area in Argentina is indicated in Figure 1. Most of Argentina's 
soybeans are grown in a concentrated area in central Argentina where the 
climate is predominantly humid subtropical. The western edge of the soybean 
area in central Cordoba is semi-arid with warm to hot summers. There, drought 
and high temperatures can be a problem during the growing season. The growing 
season begins with planting in October and November and extends to harvest from 
April through June. 


METHOD 

Multiple regression analysis of yield with selected agroclimatic indices 
was used to derive a suitable model. The regression equation for the soybean 
model derived is: 

Y = a+ BiT + B2(RDFN1) + B3(P-PET2) + E 

where 

A 

Y = Estimated yield 
cc = Constant 

Bj = Coefficients of the variables j = 1 - 3 
T = Trend 

RDFNl = Deviation from normal of the precipitation for January 
P-PET2 = Precipitation minus PET (potential evapotranspi ration) for February, 
and 

E = Unexplained error. 

Trend between the years 1969 and 1978 was determined from a plot of yields 
as shown in Figure 2. Both the weather variables chosen reflect soil moisture 
availability for plant growth. A more complete definition of P-PET can be 
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Figure 1. DENSIIT OF PLANTED AREA IN SOYBE.ANS 1977/78 CROP YEAR 
(Source: "Agronomic Characterization of the Argentina 

Indicator Region") 
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Figure 2. Plots of Production and Yield Versus Year for 
Argentina Soybeans. 
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found in the Appendix. Large positive P-PET values suggest wet conditions. 

In developing the models, various procedures of the Statistical Analysis 
System (SAS Institute, Inc., 1979) were used. The procedures used and the 
operations performed with each are summarized in the Appendix. The selected 
model had the highest and included variables that were significant at the 
10 per cent level and agronomical ly meaningful. 

DATA 

Crop data for Argentina from 1964 to 1979 were obtained from the Foreign 
Agricultural Service (Alan Vandagrith, personal communication, 1982). A data 
set was created with the year of yield as year of planting. Since the growing 
season spans two years, the crop data were expressed in terms of "year + 1" 

(or harvested year), and the meteorological data corresponding to any yield 
included data for that year and "lagged data" for the previous year. 

The meteorological data was created using the general Argentina meteorolog- 
ical station file. This file was composed of data from several different 
sources, including the Monthly Climatic Data for the World and the Servicio 
Meteorological Nacional in Argentina (R.E. Jensen, C.M. Sakamoto, and 
S.E. Mummert; August, 1974). Stations inside the soybean-producing area were 
averaged by province and weighted according to the percentage that their 
province contributed to the total country production. Figure 3 shows the 
soybean-growing area in Argentina and associated stations; Table 1 lists the 
stations used and their weights. 

The years between 1965 and 1976 (years of harvest) were used to model since 
these years contained the most complete meteorological data. 

PROCEDURES 

Weather variable selected from correlations and plots with yield. These 
variables along with TREND were then tried in regression equations. The same 
models were also tried in regression with the residuals of yield after the 
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Argentina Soybean Model. 
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STATION NAME 


WEIGHT 


BUENOS AIRES .08 

Pergamino 

Junin 

Nueve de Julio 

CORDOBA . 1 7 

Villa Dolores 
Pilar 
Bell Ville 
Rio Cuarto 

SANTA FE . 75 

Parana 
Bell Ville 
Rosario 
Casilda 
Pergamino 


Table 1. Meteorological Stations Used for 
the Argentina Soybean Model. 
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effects of trend had been removed (PROC 6LM) as the dependent variable. The 
result was that the models in which trend was included as a variable faired 
better an- nao a nigner coefficient of determination (r 2) in every case. It 
was oecioed to keep TREND as a variable. 

The limited number of data years severely limits the number of variables 
that can be Included in a regression model by limiting the number of degrees 
of freedom. This influenced the choice of model considered "best" for soybeans. 
11 ee variables, one of which was TREND, were considered the optimum with 12 
years of oata. The other two variables which best predicted soybean yields 
are RDFNl and P-PET2. The coefficients for both were positive, which indicates 
a good moisture supply is needed during flowering. Table 2 is a summary of the 
statistics of the model. 

TEST RESULTS 

A jackknife test was run on the model. In this test a year was eliminated 
from the crop data and the model was used to predict that year's yield. This 
process was done for each successive year from 1965 to 1976. The results show 
that the model predicted well as a whole. Test results are printed in Table 3 
and plotted in Figure 4. 
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APPENDIX 

De*inif’on 0*" Va'''a^1es 

oreci 31 tdtion minus potential evapotranspiration, is an index used 
to measure tne amount of moisture available for plant growth. Potential 
evapotranspiration is determined by the procedure developed by Thornthwaite 
(iy48) which uses only temperature: 



Where I » heat index, which Is the sum of the 12 monthly indices i, 



T = monthly temperature in ®C, and 

a = an empirical exponent, 6,75 x - 7.71 x lO’^I^ + 1,79 x 10"2 

+ 0.49. 

The duration of daylight is used to adjust potential evapotranspiration as a 
portion of 12 hours. 


Statistical Analysis System Procedures Used 


PROC CORR 

PROC PLOT 
PROC STEPWISE 


Computes correlation and coefficients between 
variables, including Pearson product*moment and 
weighted product -moment correlation. 

Graphs one variable against another, producing 
a printer plot. 

Provides five methods for stepwise regression. 
Stepwise is useful when selecting variables to 
be included in a regression model from a 
collection of independent variables. 


PROC STEPWISE FORWARD Begins by finding the one-variable model that 

produces the highest r 2. For each of the other 
independent variables, FORWARD calculates F- 
statistics reflecting the contribution to the 
model if the variable were to be included. 
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PRUC STEPWISE BACKWARD Begins by calculating statistics for a model 

Including all the independent variables. The 
variables are deleted from the model one by 
one until all tne remaining variables produce 
F-stat1st1cs significant at tne .10 level. 

WISE STEPWISE Tne stepwise method Is a mocificatlon of the 

forward selection technique, differing In 
that variables already in the model do not 
necessarily stay there. After a variable Is 
added (as 1n the forward selection method) the 
stepwise method looks at all the variables 
already included In the model and deletes any 
variable that does not produce an F-statlstIc 
significant at the .10 level. Only after this 
check Is made and the necessary deletions 
accomplished can another variable be added to 
the model . 

PROC STEPWISE MAXR {Maximum r 2 Improvement) Unlike the three 

techniques above, this method does not settle 
on a single method. Instead It looks for the 
“best" two-variable model, the "best" three 
variable model, and so forth. 

PROC PETM Uses latitude and mean monthly temperature to 

calculate Thornthwaite's potential evapo- 
transplratlon for each month, 

PROC ZINOEX Uses monthly PET's, precipitation, SS 

(beginning moisture in surface layer), 

AWCS (available water capacity In surface 
layer), SU (beginning moisture In the 
underlying layer), and AWCU (available 
water capacity In the underlying layer) 
to calculate Palmer's soil molstu^ 
budget, drought Index Z, ET, and ET, 





5 

> 

A 

7 


u — 
A e 

i S 

a • 

a e 




u 

A 

ff 

c 

a 

a. 


CM^e 

eeo 

• • • 

coe 


u. 7'— — 

• • ■ 

f\ — 


c 








y 







>■ 






U) 





— 


V4! 

yts. 

(/■ 

fV 


a c 

a 

<3C 

tr 

J'XC 

.r 

4 e 

< 

X.-*! 




— e 

* 

ITC 


J'X — 

> 

a c 

c 

nn 


nc— 


< e 

jr 

a'cc 


j'j'ir 

if 

> e 


ore 

Ik' 

fve.j' 


o 

2 

yif 

a 

e«v>S’ 


H- • 

< 

• • 

> 

• • • 


r ^ 

u. 

cs 


•ftrx 

j" 

u. 

S 

fV 




u. 

c. 



li. 

II 




1 








1 








1 



a 





1 



— / 

mm 




1 



u. 

u 

t' 

X ■fo 

a 

1 





u 

•tflTC 

c 

1 





■> 

rj s 

a 

X 1 





« 

"VS 

a 

e> t 1 






=^ir> 

ui 

i 



* 



<f f\) 


9xC4 1 

• 


H 


tr 

L''»crv 

c 

COC 1 

•J 


c. 



e c c 

k* 

— 9C 1 

u 


c 


U. 

• • • 


• • • 1 

> 


a 


mm 



09C 1 



a 

Km 


xC 


1 

t 







1 



*. 

r 

• 



1 

e 




7 



1 



M 





1 

c 



a 




1 


— 


• 




1 

• 



s 




1 

o 


mm 





1 



a 

M 




1 

oj 

— 





u 

r*— fr.s 1 

Z 



U 




O I 


— 





ll 

'Nj— I 




« 



«r 

nvcxC 1 

»» 

j* 


w 



> 

.rs xf* 1 

a 







1 


— 

< 




X 

r 1 

few 


y 


u 

»^a — 


r"j-C9 1 

• 

• 



A“ 

»• 


• • * • 1 

e 


U 





—• Z~ X 1 

C' 

/ 

< 





M 1 








1 

u. 







• 

mm 







1 

z 







1 

LT 







1 

s«« 







I 

4/ 



Aik 




1 




u 




I 

u 







1 




u 


c 


»- 1 

i* 



I* 


•■V 


C t 




z 


X 


U.' 1 

1 



u 


X 


U' f\. 1 






U 7 


2 r—i- 1 




C' 


r f < 


u rrii 1 

3^ 





cai- 


i^ikua 1 

2 





u t f 


r ac 1 








-*-aa 1 

u» 



< 




( 

X 



•m 




1 

i» 



0 




1 




< 




1 

X 



> 




1 

A 

b 







1 

1 

a 







1 

UJ 



< 




1 

wl 


ORIGINAL PAGE 19 
OF POOR QUALITY 


-srsT^ 


Table 2. Statistics of 'rgentioa Soybean Model. 
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Table 3. Jackknife Test Results for Argentina Soybean Model. 
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Figure 4. Argentina Soybean Model 



